Soft elastic surfaces as a platform for particle self-assembly 
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We perform numerical simulations to study self-assembly of nanoparticles mediated by an elastic 
planar surface. We show how the nontrivial elastic response to deformations of these surfaces leads 
to anisotropic interactions between the particles resulting in aggregates having different geometrical 
features. The morphology of the patterns can be controlled by the mechanical properties of the 
surface and the strength of the particle adhesion. We use simple scaling arguments to understand 
the formation of the different structures, and we show how the adhering particles can cause the 
underlying elastic substrate to wrinkle if two of its opposite edges are clamped. Finally, we discuss 
the implications of our results and suggest how elastic surfaces could be used in nanofabrication. 



I. INTRODUCTION 

Elastic surfaces are ubiquitous in nature and technol- 
ogy and appear across all length scales, from the cellular 
microcnvironment to large-scale objects such as bridges 
and buildings. The mechanical properties of these sur- 
faces play an important role both in their biological func- 
tion, and in their wide usage in material engineering. For 
instance, it is known that the stiffness of an elastic sub- 
strate alters the morphology and dynamics of tissue cells 
adhering onto it [T]. Variable cytoskeleton assembly [2] 
and cell spreading [5] on substrates of different mechani- 
cal properties are two nice examples of this. Furthermore, 
their response to external stress have been exploited in 
metrology [U |S] and in the production of micro- and 
nano-scale patterned surfaces that may serve as compo- 
nents with novel optical, electronic and magnetic prop- 
erties [6]. 

We are interested in understanding how elastic surfaces 
could be used to template aggregation of nanocompo- 
nents. The idea of using interfaces, specifically fluid ones, 
as a means of driving self-assembly of arbitrary building 
blocks was first introduced by Whitesides et. al. [HIH]. 
In this case, the presence of floating objects on a fluid 
interface induces local deformations in its profile which 
are minimized when the objects are isotropically driven 
close to each other [9l [10]. By manipulating the tension of 
the interface, and by tailoring the chemistry of the build- 
ing blocks, millimeter-size objects, microchips and micro- 
components have been successfully self-assembled [TT]- 
[13]. Unlike fiuid interfaces, whose large scale physical 
properties are dominated by their surface tension, elastic 
surfaces resist stress and respond to it in a spring-like 
fashion. 

Elastic surfaces bend and stretch in response to de- 
formations. The resulting macroscopic behavior is char- 
acterized by strong nonlinearities ^14] . The mechani- 
cal properties of macroscopic elastic sheets have recently 
been the subject of intense investigation [H irSHTO] . Un- 
der an applied force an elastic surface deforms in a way 
that minimizes the energy associated with the deforma- 
tion. It is easy to show [T31 [H] that the ratio between 



stretching and bending energies for an arbitrary deforma- 
tion of amplitude ft, on a flat elastic sheet of thickness t 
scales as Eg/Eb ^ (h/t)'^. Therefore, for sufficiently thin 
sheets, bending is the preferred mode of deformation and 
unstretchability can be thought of as an overall geomet- 
rical constraint to the deformations. The net result is 
that thin elastic surfaces respond to an external applied 
stress with stretch-free deformation involving (when pos- 
sible) exclusively uniaxial bending. Skin wrinkling under 
applied stress [TH] and stress focusing via d-cone for- 
mation of crumpled paper |19j are two examples of this 
phenomenon. 

Such nontrivial phenomenology extends to the micro- 
scale. There are several artificial and naturally occur- 
ring examples of microscopic elastic surfaces, including 
graphite-oxide sheets [201 121] , graphene sheets [211 123] , 
cross polymerized biological membranes [21], cross poly- 
merized hydrogels [3], buckypaper [25112 7j . the spectrin- 
actin network forming the cytoskeleton of Red Blood 
Cells [211 [21], and very recently they have been fabri- 
cated using close-packed nanoparticle arrays [5D]. Our 
expectation is that diffusible particles adhering over an 
elastic surface should be driven to aggregate into config- 
urations that reduce the mechanical cost of the overall 
surface deformation. These configurations will depend 
on the geometry of the surface, its elastic properties and 
the strength of the adhesion (the applied force). 

We have recently shown how the response to defor- 
mations of elastic nanotubes can alter the elastic prop- 
erties of a flexible fllament binding to it [31], and that 
elastic nanotubes and nanoshells can drive self assembly 
of nanoparticles in a variety of patterns that depend on 
the interplay between bending and stretching rigidities of 
the template, and the amount of deformation of the sur- 
face [3U 133] • In this paper we explore the phase behavior 
of nanoparticles adhering onto a planar (extended) elas- 
tic substrate as a function of the mechanical properties of 
the substrate, namely its stretching and bending rigidity, 
and the strength of the adhesion. We also analyze the 
role of the boundaries of the elastic sheet and their influ- 
ence on the aggregation patterns. Our findings suggest 
that the geometrical features of the anisotropic aggrega- 
tion of the particles can be tuned in a variety of patterns 
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by controlling the elastic parameters of the problem. 



II. METHODS 

The elastic plane is modeled via a standard triangu- 
lated mesh with hexagonal symmetry |34j. To impose 
surface self-avoidance we place hard beads at each node 
of the mesh. Any two surface beads interact via a repul- 
sive truncated-shifted Lennard- Jones potential; 
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where r is the distance between the centers of two beads, 
a is their diameter, and e = lOOfceT. 

We enforce the surface fixed connectivity by linking 
every bead on the surface to its first neighbors via a har- 
monic spring potential 



U stretching — Ks{r ~ Vb)'^ 
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Here Kg is the spring constant and it models the stretch- 
ing rigidity of the surface, r is the distance between two 
neighboring beads, rg = 1.23cr is the equilibrium bond 
length, and it is sufhciently short to prevent overlap be- 
tween any two triangles on the surface even for moderate 
values of Kg. 

The bending rigidity of the elastic surface is modeled 
by a dihedral potential between adjacent triangles on the 
mesh: 
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where (fi is the dihedral angle between opposite vertices of 
any two triangles sharing an edge and Kb is the bending 
constant. 

Particles of diameter CTc — 10a are described via the 
repulsive truncated-shifted Lennard- Jones potential in- 
troduced in Eq. [T] with cr — >■ CTc- The generic binding 
between the nanoparticles and surface is described by a 
Morse potential: 
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where r is the center-to-center distance between 
a nanoparticle and a surface-bead, r^s is bead- 
nanoparticle contact distance r^vs = 5.5a and Dq is the 
binding energy. The interaction cutoff is set to lOcr and 
7 = 1.25/cr. 

The simulations were carried out using the LAMMPS 
molecular dynamics package |35) with a Langevin dy- 
namics in the NVT ensemble. Dimensionless units are 
used throughout this paper. The timestep size was set to 
dt = 0.002to (tq is the dimensionless time) and each sim- 
ulation was run for a minimum of 5 TO® iterations. In this 
study we considered unconstrained and edge-constrained 



sheets. To minimize edge effects in unconstrained sheets 
we considered surfaces with an overall circular geome- 
try. Two different equilibrium radii Rpiane — 50a and 
Rpiane — 60.4(7 wcrc explored. To preserve the me- 
chanical stability of the sheet the colloids were placed 
both on top and at bottom of the surface. When edge- 
constrained surfaces were considered, a rectangular shape 
was selected and the colloids were placed only on one 
side of the plane. For this specific case we considered 
two surface equilibrium areas, A ~ (176 x 152)a^ and 
A ~ (244 X 212)cr^. In both cases a wide range of 
nanoparticle surface fractions between 10% and 60% was 
explored. Typical values of ct ~ 10 — 20nm would imply 
colloidal particles of diameter ^ 100-200 nm and surfaces 
of area A ^ 200-1000/xm^ Figure[l]illustrates the model 
used in our simulations. 



III. RESULTS 

We find that elastic surfaces can drive nanoparticle ag- 
gregation. The geometry of the aggregates can be tuned 
into a variety of patterns controlled by the mechanical 
properties of the surface {Kg and Kb) and the strength 
of the particle's adhesion (Dq). Let us begin by looking 
at the role of the membrane's stretching rigidity. Fig.[T^) 
shows a diagram of the different aggregates obtained for 
different values of Kg as a function of the extent of the 
surface deformation (regulated by Dq) at fixed bending 
rigidity Kb- Fig. [l|;) shows simulation snapshots of the 
corresponding patterns. As expected, when Dq is small 
the surface is basically unaffected by the presence of the 
particles and the particles behave effectively as a low- 
density two-dimensional hard sphere fluid. In the op- 
posite limit, when the particles bind very strongly, the 
membrane undergoes large local deformations limiting 
the diffusion of the particles and resulting in kinetically 
trapped configurations. Repeating the simulations under 
the same conditions leads to a different not well defined 
configuration. We call this phase the arrested phase. 

The intermediate regime is characterized by five dis- 
tinct structured phases. For small values of Kg the ag- 
gregation is completely driven by the minimization of 
the bending energy. As a result particles aggregate into a 
two-dimensional hexagonal crystal. This is what happens 
for instance in lipid bilayers where a bending-mediatcd 
isotropic, V ^ —r^"^, interaction can drive surface inclu- 
sions at a separation r from each other to aggregate [36] . 
Upon a small increase in Kg the hexagonal crystal rear- 
ranges into a square lattice. We believe the coexistence of 
two lattices is either due the substrate having the elas- 
tic modulus near the boundary between the hexagonal 
and square crystalline phases, or due to small spatial in- 
homogeneity of the substrate's elastic modulus, either 
way, this recent experimental observation seem compati- 
ble with our numerical results. As Kg is further increased 
we find that the crystalline aggregate is disrupted in favor 
of a network of short connected lines. This phase is the 



3 



consequence of balancing the stretching and the bending 
energy: the former preferring the stretch-free uniaxial 
deformations, and the latter driving in-plane isotropic 
aggregation. 

For even larger values of Ks the connected network 
is disrupted and particles arrange into straight parallel 
lines. Increasing Ks at this points only leads to to a 
larger stiffness of the linear aggregates. This transition is 
completely stretching-driven. The parallel lines start ap- 
pearing when Kg ^ K^- This is clearly shown in Fig. [T|d) 
where we show how the formation of straight and con- 
nected aggregates depend on both stretching and bending 
constants. For Kg/Kb ^ 1 one indeed recovers the thin 
and unstreachable sheet limit for which only stretch-free 
(uniaxial) deformations are possible. The most dramatic 
consequence of this property of elastic plates is the fifth, 
folded phase. This phase occurs for larger values of Dq, 
when particles tend to increase the contact area with 
the membrane as much as possible. In this region the 
surface immediately folds into a well organized higher 
three-dimensional hexagonal structure (Fig.jlJ:)). 

To better characterize the dependence of the different 
phases on Kg - from the hexagonal to the square lattice, 
from the connected network to the linear one, we also 
measured the frequency of particle contacts as a function 
of Kg. Fig. [2] shows the probability distribution of the 
number of the nearest neighbors as a function of Kg in 
the different phases. The connectivity decreases when 
increasing Kg, going from the six neighbors of the hexag- 
onal phase to the four neighbors in the square lattice, 
and finally two for the connected and the straight lines. 
For the linear aggregates the significant difference is not 
in the location of the peak of the distribution (indeed a 
large number of particles will have two neighbors even 
in the connected linear aggregates), but in the relative 
height of -P(3) which is the signature for branching. 

It should be emphasized that the number of connec- 
tions does decrease continuously with increasing Kg. It 
is tempting to interpret these data in terms of a sin- 
gle growing length scale that sets the size for the aver- 
age distance between any two nodes in the linear net- 
work, and consider the straight-line phase as the limit- 
ing behavior in which this distance becomes larger than 
the surface. A simple mean field calculation balancing 
stretching and bending energies [19| points to the length 
scale Ip cx h^/'^{Kg/ K^^Y/'^R^/'^ , which qualitatively pro- 
duces the correct phenomenological behavior, but unfor- 
tunately the small system sizes analyzed in this study 
prevent us from making such a link more concrete. 

It is important to stress that the free boundaries of the 
membrane play an important role. Indeed, it is not clear 
whether the linear phases indicated in Fig. [l] are stable 
with respect to folding. In fact, in a few cases our longest 
simulations of the linear phases resulted eventually in 
a folded phase. We expect this to be an effect of the 
free boundary of the surface that can be taken care of 
by applying a small external tension or by clamping the 
outer edge of the surface. To show that this is indeed the 



case, we considered a rectangular elastic sheet in which 
two opposite sides (edges) are kept fixed (clamped). In 
the absence of the adhering particles the sheet remains 
flat to its equilibrium size. Once the particles bind to it 
we observe only two phases for moderate values of Kg, 
the gas phase and the straight-linear phase. The former 
appears when Dq is insufficient for the particles to deform 
the membrane, while the latter occurs when crosses 
a certain threshold value which depends mainly on the 
bending rigidity of the plane (Fig.|3^)). 

The linear structures formed in this phase always ap- 
pear to be perpendicular to the constrained sides of the 
membrane (Fig. |3]3)). However we observe that the dis- 
tance between them can be tuned by changing Kg and 
Ki,. This kind of pattern is reminiscent of the wrinkle 
pattern that occurs when a thin elastic sheet is subjected 
to a longitudinal stretching strain[T6l [17]. The sheet 
is then unable to contract laterally near the clamped 
boundaries, so it wrinkles to accommodate the in plane 
stress. Cerda and Mahadevan showed that, for a con- 
stant tension, the wavelength of the wrinkles scales as 
A ~ (L)i/2 ^ (if,/X,)i/4 [T7]. 

Indeed, we find the same reasoning can be applied here. 
Instead of having an external force stretching the plane, 
the particles binding to the surface act as the stress 
source causing the sheet to wrinkle. This stress is di- 
rected perpendicularly to the fixed sides of the plane. 
Since the wrinkles are the regions where the particles 
can gain the highest contact area with the surface, i.e the 
highest binding, the particles follow the wrinkle pattern 
resulting in the straight parallel aggregates. We believe 
that the destabilization of the linear-connected phase is 
due to the implicit symmetry breaking imposed by the 
way we clamp the membrane, in fact, when clamping is 
enforced on all four edges of the sheet, the phase reap- 
pears. 

We also analyzed the dependence of the wavelength of 
the particles' lines with Ki, and Kg and it appears to 
nicely follow the theoretical prediction of Cerda and Ma- 
hadevan (Fig. ^p))- Nevertheless, two extra parameters 
play a role in determining the separation between the 
lines in this case: the surface coverage and the particles's 
binding energy. Since particle binding to the surface is 
favorable, once the particle's density becomes larger than 
that required to completely fill the wrinkles with parti- 
cles, new lines (wrinkles) form in between the preexisting 
ones, bringing the preexisting ones closer together. The 
inset in Fig. [sja) shows the decrease in A with the in- 
crease in the particle density, for two different values of 
Kg. In addition to that, we find that the increase in Dq 
(for constant Kg and K^) also brings the lines closer to- 
gether. Higher binding increases the amplitude of the 
wrinkles (analogous to increasing the strain tension in 
[IB]), which decreases the surface area accessible to the 
particles, effectively increasing the density. 

It should be stressed that the mechanism driving 
self-assembly of particles into linear aggregates that we 
described is significantly different from the controlled 
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wrinkling methods recently developed for the fabrica- 
tion of patterned surfaces [6]. There the wrinkles are 
preformed by compressing the substrate, and particles 
trivially arrange along the wrinkles' axis to maximize 
their binding energy, in our case the surface is not 
pre-wrinkled, and the linear aggregates develop (in 
a reversible manner) as a result of a more delicate 
balance between the energies of the system and the 
collective behavior of the particles. Interestingly, once 
the wrinkled phase is formed it is possible to control the 
overall direction of the lines by simply applying a small 
external tension. For instance, the release of the surface 
clamping and simultaneous application of a small tension 
in the direction perpendicular to the direction of the 
wrinkles, results in a reorientation of the lines along the 
direction of the tension. This supports our assumption 
of the reversibility of the line-forming process and 
could suggests even richer potential application of this 
approach for periodical patterning. 



CONCLUSIONS 

In this paper we shown how elastic surfaces can tem- 
plate self-assembly of nanoparticles, similarly to the way 
fluid interfaces do. We show how by tuning the relative 
cost of bending and stretching energies (i.e. the thick- 
ness of the sheet) it is possible to control the geometry 
of the aggregates. The formation of the different lin- 



ear aggregates, for thin sheets, is an explicit manifesta- 
tion of the anisotropic interaction between the nanoparti- 
cles. When the surfaces become effectively unstreachable 
particles arrange into macroscopic ordered parallel lines 
whose separation can be controlled by the elastic param- 
eters of the surface. Clamping of the edges across the 
membrane substantially improves the periodic ordering 
in the system. 

The physical properties of our model can be mapped 
onto a model of a thin sheet supported on an elastic foun- 
dation if the stretching rigidity of the plane is substi- 
tuted by the stiffness of the elastic foundation. There- 
fore, the results of our theoretical study are quite gen- 
eral and may suggest novel use of the elastic interfaces 
in nano/micromechanics and material engineering. Pos- 
sible experimental systems where our predictions could 
be tested include cross-polymerized or crystalline lipid 
bilayers, thin polymeric sheets, ultrathin cross-linked 
nanoparticle-membranes or possibly free standing liquid 
crystalline films in the presence of colloidal particles. 
More in general, any elastic substrate that can be lo- 
cally deformed by the interaction with a diffusable bind- 
ing component. 
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FIG. 1: Illustration of the triangulated mesh model used in our simulations. The surface beads of diameter a (blue spheres) 
are set at the nodes of each triangular element to enforce surface-self-avoidance and are linked to their first neighbors with 
springs of the constant Kg and the equilibrium length rjvs (measured form the beads centres). The surface connectivity is kept 
constant, and apart from boundary nodes each surface bead has six neighbors. The dihedral angle 1-2-3-4 from which bending 
energies are computed is also indicated. This energy is minimized when all angles between neighboring triangles are equal to 
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FIG. 2: (a) Phase diagram of nanoparticles binding to an elastic planar surface. In this case the bending rigidity is Kb = ISOfcaT; 
the equilibrium radius of the surface is i? ~ 60.4(t and the number of nanoparticles is A'' = 40. The lines separating the different 
phases serve as a guide to the eye. The arrow points in the direction of lower line connectivity, (b) Boundary between connected 
to straight parallel lines as a function of Ks and Kb- The dashed Kb — Ks line serves as a guide to the eye. (c) Simulation 
snapshots of the seven observed phases. For the sake of clarity the hexagonal and the square crystal phases are shown with the 
larger number of particles than the other phases. 
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FIG. 3: Probability distributions of particle contacts in a self-assembled aggregates for different values of Ks, and constant 
Kb = 150fcBr. Prom right to left the distributions refer to the planar hexagonal crystal, the planar square crystal, the 
interconnected lines, and the straight parallel lines. 
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FIG. 4: (a) Phase diagram of nanoparticles binding to a clamped rectangular elastic surface. Two phases are observed for 
different values of Kt and depend only on Dq: the gas phase and straight parallel lines phase (SL). These data refer to 
the case in which Ks = 150kBT/a^, N = 60, and the area of the the plane is (176 x 152)cr^. (b) Line separation A as a 
function of the mechanical properties of the surface. We show data for two different surface coverages: (j) = 38.7% (cross 
symbols) and </!> = 21.2% (circle symbols). The straight lines represent the fit of the data to the scaling law A ~ (Kb/ Kg)^^^ . 
The inset shows the dependence of A on the particle surface coverage, shown for three different values of the {Kb/ KsY^'^ 
parameter: 1.35CT^^^(top), 0.9cr^''^ (middle), O.SIct^''^ (bottom), (c) Simulation snapshots of linear aggregates for three different 
combinations of the elastic parameters and the surface coverage densities. The equilibrium surface area is A = (176 x 1^2) . 



